Frustrated spin ladder with alternating spin-1 and spin- 1/2 rungs 
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We study the impact of the diagonal frustrating couplings on the quantum phase diagram of a two- 
leg ladder composed of alternating spin-1 and spin-1/2 rungs. As the coupling strength is increased 
the system successively exhibits two gapped paramagnetic phases (a rung-singlet and a Haldane- 
like non-degenerate states) and two ferrimagnetic phases with different ferromagnetic moments per 
rung. The first two states are similar to the phases studied in the frustrated spin-1/2 ladder, whereas 
the magnetic phases appear as a result of the mixed-spin structure of the model. A detailed char- 
acterization of these phases is presented using density-matrix renormalization-group calculations, 
exact diagonalizations of periodic clusters, and an effective Hamiltonian approach inspired by the 
analysis of numerical data. The present theoretical study was motivated by the recent synthesis 
of the quasi-one-dimensional ferrimagnetic material Fe 7I Fe 7I/ (trans- 1,4-cyclohexanedicarboxylate) 
exhibiting a similar ladder structure. 

PACS numbers: 75.10.Jm, 75.50.Gg, 64.70.Tg 



I. INTRODUCTION 

Over the past two decades there has been an increas- 
ing interest in quantum spin systems with competing 
exchange interactions^^ Quantum spin chains and lad- 
ders with frustration, both for half-integer and integer 
spins, set up an important part of this research since 
they provide a unique testing ground based on the avail- 
able powerful analytical and numerical techniques for 
one-dimensional (ID) systems. In particular, the frus- 
trated ladder models have allowed controlled calculations 
to examine topological order^ dimer order as well as 
the appearance of fractional excitations in spin models^ 
Most of previously studied frustrated chain and ladder 
models have been related to uniform-spin structures with 
all the spins same. In comparison, till now much less 
experimental as well as theoretical work concerning the 
impact of competing interactions in quasi-lD mixed-spin 
systems has been accomplished £ Often these systems ex- 
hibit quasi-lD ferrimagnetic ground states with a net 
ferromagnetic moment, so that apart from rich quan- 
tum phase diagrams they might be expected to provide 
generic examples of ID magnetic-paramagnetic quantum 
phase transitions^ 

On the experimental side, during the past two decades 
it has become possible to synthesize a large variety of 
quasi-lD materials with ferrimagnetic properties. Most 
of these materials are heterometallic molecular magnets 
containing different transition metal ions in the unit celli^ 
A generic spin model describing these materials is the 
quantum Heisenberg spin chain with antiferromagnetic 
nearest-neighbor exchange interactions and two types of 
alternating quantum spins with magnitudes Si and S2 
(Si > S'2)^— In the extreme quantum case of spins 
(1, 1/2), the latter model was shown to provide an ex- 



cellent description of the thermodynamic parameters of 
the recently synthesized quasi-lD bimetallic compound 
NiCu(pba)(D 2 0) 3 -2D 2 (pba = 1,3-propylenebis)^ An- 
other important class of quasi-lD ferrimagnets - the so- 
called topological ferrimagnets - is related to some ho- 
mometallic materials exhibiting composite chain struc- 
tures with different magnetic sublatticesJ^ The ho- 
mometallic material A3Cus(P04)4 (A=Ca,Sr,Pb) is an 
example of such quasi-lD ferrimagnets: In this com- 
pound, the Cu 2+ ions form diamond chains with strongly 
coupled trimers bridged by oxygen ions^ Since quasi-lD 
homometallic materials usually have rich exchange path- 
way structures, they may be expected to provide some 
real examples of quasi-lD ferrimagnets with magnetic 
frustration. To the best of our knowledge, the recently 
synthesized mixcd-valcnt magnetic material Fc H Fe HI 
(trans-l,4-cyclohcxancdicarboxylate)i£ provides the first 
real example of a quasi-lD Heisenberg ferrimagnet with 
magnetic frustration^ The experimentally established 
magnetic structure for temperatures larger than 36 K 
corresponds to the mixed-spin ladder with diagonal ex- 
change bonds shown in Fig. [H where the site spins 
Si = 5/2 and S2 = 2 are respectively related to the 
magnetic ions Fe //J and Fe 11 

The mentioned experimental achievements motivated 
a series of theoretical studies on quantum mixed-spin 
chains and ladders with geometric frustration. The 
symmetric diamond chain with antiferromagnetic ver- 
tical bonds was probably the first studied model of a 
ID quantum ferrimagnet with competing interactions^ 
A variant of this model, the distorted spin-1/2 dia- 
mond chain, has received special theoretical^ as well 
as experimental^ interest due to its rich quantum 
phase diagram^ and the relevance for the real mate- 
rial Cu3(C03)2(OH)2. The diamond Heisenberg chain 
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FIG. 1: The mixed-spin ladder considered in the paper. The 
arrows show the classical canted state described by the an- 
gles < <j> < tv/2 and < 8 < 7r/2 for the classical spins with 
magnitudes Si and 52, respectively. The other two classical 
phases correspond to spin configurations with (<j), 8) = (0, 0) 
(antiferromagnetic state) and (<j>,8) = (71-/2,71-/2) (ferrimag- 
netic state). 

is also one of the simplest quantum spin models admit- 
ting four-spin cyclic exchange interactions^ A generic 
quantum spin model of a frustrated ID ferrimagnet is 
the mixed-spin Heiscnbcrg chain composed of two types 
of alternating spins interacting via competing nearest- 
neighbor and next-nearest-neighbor antiferromagnetic 
exchange bonds*2I This model may also be considered 
as a mixed-spin zigzag ladder and is a ferrimagnetic 
analogue of the frustrated Heisenberg chain with ferro- 
magnetic nearest-neighbor and antiferromagnetic next- 
nearest-neighbor exchange bonds. The spin-1/2 frus- 
trated Ji — J 2 ferromagnetic chain has recently at- 
tracted much attention^ as it is supposed to describe 
a number of quasi-lD edge-sharing cupratcs, such as 
Rb 2 Cu 2 Mo 3 Oi2^, Li 2 ZrCu0 4 r^ and LiCuV0 4 ^ The 
latter material exhibits multiferroic properties^ as well 
as an interesting specific phase transition in a magnetic 
field from an ordered spiral to an ordered modulatcd- 
collinear magnetic phases j^I There are other two generic 
types of frustrated mixed-spin ladder models describing 
two interacting mixed-spin alternating chains. The first 
one is the checkerboard mixed-spin Heisenberg ladder 
with frustrating diagonal exchange couplings^ and the 
second one is the two-leg ladder model with two types of 
alternating rungs presented in Fig. Q] Finally, there has 
been a lot of recent work reporting interesting quantum 
phase diagrams in different composite Heisenberg chains 
with ferrimagnetic ground states^ 

In this study we focus on the effects of frustration on 
the ground state phase diagram of the mixed-spin ladder 
shown in Fig. [TJ In addition to the theoretically inter- 
esting question of the effects of frustration in this sys- 
tem, an experimental realization of a closely related sys- 
tem in a mixed-valence iron polymer further motivates 
usJ£ In the next section we introduce the model and 
study some relevant properties of its Hamiltonian. In 
Section [TTI1 we give a detailed description of the quan- 
tum phases by using an effective Hamiltonian approach 
inspired by the analysis of data obtained using density- 
matrix renormalization-group (DMRG) and exact diag- 



onalization (ED) techniques. We conclude in Section HVl 
with a brief summary of the results. 

II. THE MODEL 

The system under consideration (sec Fig. [T|) consists of 
two equivalent mixed-spin Heisenberg chains (character- 
ized by the nearest-neighbor exchange constant J3 > 0) 
coupled via rung ( J\, J 2 > 0) as well as diagonal (J4 > 0) 
exchange bonds. The Hamiltonian of the system reads 
as 

T~L = H12 + T~iz + T~ii, (1) 

where 

L/2 

T~i\2 = ( Jl<Sl,2n ' S2,2« + J 2 <T l,2n-l ' tT 2 , 2? t-l) , 

L/2 2 

H3 = J3 /J [s m ,2n ' (c?ri, 2n-l + &m,2n+l)] , 

n— 1 m— 1 
L/2 

Hi = Jj [si,2n " (o"2,2n-l + P"2,2n-n) 
(<Tl,2n-l + <?"l,2n+l)] • 

Here Sfc j2rt and <Tfc, 2n -i (k = 1,2) are, respectively, spin- 
s'! and spin-S^ operators (Si > <S 2 ), and L is the number 
of rungs. 

It is instructive to present the Hamiltonian in the fol- 
lowing form 

L/2 

s 2n ' ( cr 2n-l + cr 2n+l)] + J aV , (2) 

n=l 

where J s , a = (J 3 ± Ji)/2, and s 2n = «i,2n + «2,2n and 
<T2n+i = 0%2n+i +""2,2n+i are rung spin operators. The 
operator V reads as 

L/2 

V = = I/ 2 „ • (hn-l + hn+l) , (3) 

where L 2n = Si,2n - s 2,2n and Z 2 n±i = Ci,2n±i - o" 2 ,2n±i 
are rung vector operators. The following analysis of the 
zero-temperature quantum phase diagram addresses the 
extreme quantum case of spins Si = 1 and S2 = 1/2, and 
is mainly restricted to the parameter subspace defined by 
J 1 = J 2 = J 3 > and J4 > 0. To some extent, such a 
choice of the parameters is motivated by the experimen- 
tally established strengths of the exchange couplings in 
the ferrimagnetic ladder material Fc 77 Fe 777 (trans-1,4- 
cyclohexanedicarboxylate) .— 

A. Symmetries of the model 

The mixed-spin system inherits some important sym- 
metries of the parent uniform-spin Heisenberg ladder 
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with diagonal interactions^ First, if the parameters 
J 3 and J 4 in % are exchanged, one can recover the 
original Hamiltonian by exchanging cither the spins on 
the Si rungs (si,2n < — > s 2,2n), or the spins on the 
52 rungs (or^n-i < — > f2.2n-i)- This means that 
H(Ji, Ji, J3, Ji) = H( Ji, J2, Ji, ^3)- Therefore, the 
study of the model can be restricted in the region J4/J3 < 
1 since the model with </ 4 / J3 > 1 maps onto the one with 
Ji/ 'h < 1- Because of the same symmetry, the Hamilto- 
nian @ does not contain mixed products of rung spins 
and rung vector operators. 

The second property of H concerns the subspace J3 = 
Ji {J a = 0), when the last term in Eq. fl2J) disappears. 
As is the uniform-spin case,— in this parameter subspace 
the Hamiltonian H commutes with the local operators 
s 2n and <y\ n _ x (n = 1, 2, . . . ,L/2), which means that the 
rung spins s 2n and a 2n -i [defined as s\ n = s 2n (s 2 „ + 1) 
and 0"2n-i = o"2re-i(CT2n— 1 + 1)] are good local quantum 
numbers. Thus in every sector of the Hilbert space, de- 
fined by the sequence [a\, 82, ■ ■ ■ , sl), the first two 
terms in Eq. ([2]) reduce to the constant 

E Q = ~[JiSi(Si + l) + J a S 2 (S2 + l)] 

L/2 

+ 2^2 [^l S 2n(s 2n + 1) + J202„-l(o-2 n -l + 1)] • 

n=l 

Thus Eq. @ takes the simple form of a Heisenberg spin 
chain 

L/2 

H = Eo + ^Js S 2n ■ (<T2n-l + ^2n+l) ■ (4) 

n=l 

The above expression for Eq implies that for strong 
enough rung interactions (J1/J3, J2/J3 ^ 1) the sin- 
glet eigenstate of Eq. Q, defined as a product of local 
rung-singlet states, becomes an exact ground state of the 
model. This state belongs to the sector [0,0,..., 0, 0] and 
can be considered as a prototype of the rung-singlet phase 
of Eq. ^ discussed below. The following analysis of the 
quantum phase diagram of Eq. implies that in the 
extreme quantum limit (Si,S 2 ) = (1,1/2) the sectors 
[1,1,. ..,1,1], [1,2,. ..,1,2], and [1, 1, 1, 2, . . . , 1, 1, 1, 2] 
also play an important role: In the first sector, the model 
defined by Eq. ^ is equivalent to the spin-1 Haldane 
chain, whereas in the last two sectors Eq. (|4]) represents 
spin-alternating ferrimagnctic chains. The ground states 
related to these models appear in the quantum phase di- 
agram of the discussed system. 



B. Classical phase diagram 

The classical phases of Eq. ([1]) can be described by the 
angles (j> and 9 (see Fig. [TJ which determine the orienta- 
tions of the classical spins in the xz plane. We consider 
the parameter subspace defined by J\ = J 2 = J3 = 1 and 



J a > 0. The expression for the ground-state energy per 
cell containing two rungs is seen to be 

= -|l C os(20)-^cos(20) 

-4cos(0-e»)+4J 4 cos(0 + 6»). (5) 

A minimization using the independent angle variables 4> 
and 9 gives the following equations: 

cos(</> + #) = — J4 — c 2 

K 

cos(</> — 6) = C2J4 — c\K, (6) 

where c\ = a — er -1 , c 2 = o + er -1 , and a = S\jS 2 > 1. 
The parameter n = k{Ja) reads k = (4J|/3 — 1/3) 1 / 2 . 

The lower and the upper (J4 ) phase bound- 

aries of the classical canted phase shown in Fig. [1] are 
related to the inequalities | cos{4> + 6)\, \ cos(</> — &)\ < 1 
implying 
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FIG. 2: (a) The classical phase diagram described by the 
angles cj> and 9 vs. J4, as obtained from Eq. ([6]) for the system 
with Si = 1 and S2 = 1/2. (b) z components of the classical 
magnetizations in the Si (Mi) and S2 (M2) sites of the same 
system. The filled circles on the J4 axis correspond to the 
classical transition points = 7/13 and J4 ' = 1. 

For J 4 < J 4 rf \ we get states of zero magnetization 
in which the two spins on any rung and spins along a 
leg are antiferromagnetically aligned. The canted state 
realized for j| d ' < J 4 < has a net magnetization 
that takes a maximal value at some intermediate J 4 be- 
tween both boundaries [see Fig. [Hb)]. For J 4 > 
this classical canted phase gives way to a ferrimagnctic 
state where all the spins of the same magnitude are fer- 
romagnctically aligned but the relative alignment of Si 
and £2 is antiferromagnetic. Notice that the magnetic 
measurements in Ref. |l5| indicate the discussed ferrimag- 
netic configuration-eventually with a small canting of the 
classical spins- as the most probable spin configuration 
realized in the real material Fe II Fc 111 . For Si = 1 and 
S2 = 1/2, the above equations give = 7/13 w 0.538 
and J4 ' = 1. For the real material studied in Ref. [ID 
(Si = 5/2, S 2 = 2), one has J { 4 d) = 61/121 « 0.504 and 
J ( 4 u) = 21/39 w 0.553. 
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FIG. 3: Unit-cell isotropic spin-spin correlations as a func- 
tion of the frustration parameter J4, as obtained from the 
DMRG method for open boundary conditions (L = 100). 
J4 1 = 0.710, jf = 0.875, and jf = 0.975 are the special 
points identified as phase-transition points between different 
ground states. The inset shows the difference in the spin-1 
rung correlations in two neighboring cells. Note that the pre- 
sented spin-spin correlations belong to unit cells far from the 
ends. 



FIG. 4: Ground-state energy per rung as a function of the 
frustration parameter J 4 (DMRG, L = 90). JF = 0.723 
denotes the location of the maximum. The positions of the 
special points identified in Fig. [3] separate different ground 
states: RS (rung-singlet), HL (Haldane-like ), and two dif- 
ferent ferrimagnetic states (Fi and F%). The straight line ab 
represents the energy of the Haldane state (\^h)) defined by 
E H = <* ff |W|*H>. ' 



Interestingly, the discussed classical ferrimagnetic state 
appears only for relatively small values of a. For larger 
<r, the lowest energy collinear configuration for large J4 
is a non-magnetic state with ferromagnctically arranged 
legs pointing in opposite directions (i.e., antifcrromag- 
netically aligned rungs). Comparing the energies of both 

configurations (eP = Sf + 5| - 4S*iS , 2 - <iS 1 S 2 J i , 
EP = -Sf - Sf + 4S 1 S* 2 - 4S*iS' 2 J 4 , respectively), we 
see that the ferrimagnetic configuration is realized only 
in the interval 1 < a < 2 + VS « 3.73. In the large a 
case, the canted phase is also modified: On increasing the 
parameter J4 from J^ up to J4 , the S 2 spins smoothly 
change their orientation by 7r, whereas the net orientation 
of the larger Si spins coincides at the phase boundaries. 
In both variants of the classical phase diagram the phase 
boundaries are defined by Eq. (JT])- 

Finally, the discussed classical phase diagrams were in- 
dependently confirmed by our classical Monte-Carlo sim- 
ulations. Below we argue that the classical ferrimagnetic 
phase survives quantum fluctuations, whereas both the 
antiferromagnetic as well as the canted classical phases 
are completely destroyed. 



III. QUANTUM PHASE DIAGRAM 

We consider the parameter subspace defined by J\ = 
J 2 = J 3 = 1 and < J 4 < 1.5, and use the DMRG 
method^ 2 - for open boundary conditions supplemented by 
ED data for periodic clusters containing up to L = 14 



rungs. DMRG is carried out for this system for a range 
of lattice sizes up to L = 100 rungs with the spin values 
Si = 1 and S2 = 1/2, respectively. Up to 320 density ma- 
trix eigenvectors were retained. Depending on the value 
of J4, the truncation errors are between 10~ 7 and 10 -12 . 

The DMRG results presented in Fig. [3] reveal three 
special points on the J4 axis separating regions with 
different characteristics of the short-range correlations: 
J 4 cl = 0.710, Jf = 0.875, and Jf = 0.975. The same 
points are also presented in Fig. [4] which shows DMRG 
results (L = 90) for the ground-state energy of the mixed- 
spin model ([1]). A detailed numerical analysis, using both 
the DMRG and ED methods, predicts singlet ground 
states in the entire region < J4 < J 4 2 . For J4 > J| 2 , 
the same analysis suggests ground states characterized by 
net ferromagnetic moments. Below we argue that these 
special points are related to quantum phase transitions 
between different ground states. 

A. Mapping onto the frustrated spin-1/2 ladder 

An inspection of the short-range correlations presented 
in Fig. [3] implies that the weight of the local rung quintet 
(i.e., S2n = 2) states on the spin-1 rungs is negligible 
almost in the whole interval < J4 < Jf 2 . Indeed, by 
using the identity (s lt2n - s 2 .2n) = ((«L) ~ 3 / 2 ) / 2 - 5 / 4 > 
one finds that the following relation between the average 
rung correlations should be satisfied for any state with a 
zero weight of the rung quintet states: 

5 

(Sl,2n ' •S2,2n) = (f l,2n-l ' 0"2,2n-l) — T- (8) 
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As seen from the numerical results, the above relation is 
almost perfectly fulfilled in the entire region < J4 < 
J| 2 , excluding some narrow vicinity of the point J| 2 
where the correlations (si^n • S2,2n) abruptly change to 

1. The extremely small contribution of the quintet 
rung states in the region < J4 < J| 2 can be explained 
by the peculiarities of the energy spectrum of the mixed- 
spin plaquette, where the lowest quintet state happens to 
be well separated from the low-lying triplet and singlet 
states. Note that the excitation of local quintet states is 
controlled by the last term (V) in the Hamiltonian ©. 
Thus, starting from an eigenstate belonging to the sector 
S2nj o"2re-i = 0, 1 (n = 1, . . . , L/2), the first-order correc- 
tions to the wave function of this eigenstate will contain 
relatively small amount of configurations belonging to the 
sectors with local quintet states due to the larger energy 
denominator in the perturbation expression. 

These observations suggest, in particular, that in 
the discussed region the ground-state properties of the 
mixed-spin system may be approximately interpreted by 
projecting out the local quintet states in the mixed-spin 
Hamiltonian @. Up to first order in J a , the projected 
Hamiltonian reads as (see the Appendix) 



eff - —q JLi i- ^ [jj_Q"l,n • Q"2,n 

n— 1 



= —gJL 

J s cr n ■ a n+ i + 



(9) 



where o"i n and <T2 jra are spin- 1/2 operators, er n = a\ in + 

&2,n, In = Vl,n — &2,n, J± = J±, J s = Js, and J a = 

—2-^/2/3 J a . For simplicity, we have restricted ourselves 
to the case of equal rung couplings (Jj = J 2 = Jj_). The 
effective Hamiltonian © describes a frustrated spin- 1/2 
Heisenberg ladder characterized by three parameters, i.e., 
the strength of the rung (Jj_), leg (J 3 = J s + J a ), and 
diagonal ( J 4 = J s — J a ) exchange bonds. Using the same 
reasoning, it may be safely suggested that the next-order 
corrections in J a do not change substantially the singlet 
ground states, so that the effective Hamiltonian ([9]) may 
be used (i) to identify the singlet ground states of the 
original Hamiltonian © in the region < J4 < J| 2 and 
(ii) to analyze the related quantum phase transitions. 

As is well-known, as a function of the frustration pa- 
rameter J 4 the model fl5J) exhibits the so-called rung- 
singlet (RS) and Haldane-like (HL) phases ^i 3 ^ Both 
ground states are non-degenerate and exhibit finite 
singlet-triplet gaps. The character of the quantum RS- 
HL transition in the weak-coupling limit is still under de- 
bate: Some of the cited work a 30 ' 33 i 35 ' 36 suggest a direct 
first-order transition between these phases, but the oth- 
ers predict an intermediate columnar dimer phasei 4 ' 34 i 37 
Thus the mapping of Eq. ^ implies that the special 
point J4 = J! 1 can presumably be identified as a quan- 
tum phase transition point separating similar phases. Of 
course, such an analysis does not exclude the presence 
of some intermediate singlet phases in a tiny interval be- 
tween the RS and HL states. Some hints in this direction 



inspired by the DMRG results for the ground-state en- 
ergy (Fig. 2]) will be discussed below in more detail. 

The established connection with the frustrated spin- 
1/2 ladder model is additionally supported by the fact 
that the special point J4 1 perfectly maps on the RS-HL 
phase boundary in the phase diagram of the frustrated 
spin-1/2 ladder models Indeed, taking the parameters 
2/1 = Jj_/ J 3 and 1/2 = J4/J3 used in Ref. HO, the estab- 
lished relations J s = J s and J a = — 2i/2/3 J a between 
the parameters of the original and the projected Hamil- 
tonians take the form 



yi 



J±/J 3 



2/2 



6 2 - h J4/J3 
b 2 Ji/J3 ~ b\ 



(10) 



where 61 = \T2/3- 1/2 and b 2 = y/2/3 + 1/2. Note that 
the change of J4 (at fixed J± = J3 = 1) corresponds to a 
run in the (2/1 , 2/2) plane on the ab line ( see Fig. [5]) defined 
by 2/2 = (&1/&2 + 1)2/1 — &i/&2- Following Ref. H3, we 
may identify the position of the quantum phase transition 
with the point J4 = J^ 1 = 0.710 where the spin-1/2 rung 
correlations change their sign (see Fig. [3]). We find that 
the (2/1,2/2) image A of the transition point J^ 1 maps 
perfectly on the phase boundary in the (2/1, 2/2) plane. In 
Figure El we also show the symmetric point A obtained 
by the coordinate transformations y\ — > 2/1 /yi and 2/2 — > 
I/2/2, which are related to the exchange symmetry J3 < — > 
J4 of the Hamiltonian. As expected, the symmetric point 
A also lies on the phase boundary. 




FIG. 5: Phase diagram of the effective spin-1/2 ladder 
model with diagonal bonds.— The point A with coordinates 
(2/1,2/2) = (1.618,1.766) is the image of the special point 
Jf = 0.710 obtained by using Eq. pO]l. The point A is 
an image of A corresponding to the symmetry transformation 
J 3 < — > J 4 . ab is the path in the (2/1, 2/2) plane corresponding 
to the change of J4 at fixed Jx = J3 = 1. 
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FIG. 6: Unit-cell isotropic spin-spin correlations of the model 
as a function of J3 (J4 = 0, L = 100). The notations are 
defined in Fig. [3] The inset on the left shows the ground-state 
energy per rung vs. J3 and the inset on the right shows the 
variation of the singlet-triplet excitation gap with J3. 



B. Rung-singlet and Haldane-like phases 

1. Rung-singlet phase 

The RS phase, originally studied in the two-leg 
spin-1/2 ladder without diagonal bonds; 38 ' 39 is a non- 
degenerate singlet state with a finite singlet-triplet gap. 
The existence of a spin gap in this model can be eas- 
ily anticipated by using a strong-coupling analysis^ For 
J 3 / 'J_l <C 1, the ground state is a simple product of rung 
singlet bonds. The lowest rung excited states are local 
triplets with a characteristic gap cx Jj_ which survives 
the perturbation in J 3 /J ± . On the other hand, the per- 
turbation produces an energy band (with a bandwidth 
cx J 3 ) of triplet excitations. 

The same physics can be easily extracted from a 
strong-coupling analysis of the mixed-spin ladder (J2}. In- 
stead of doing this, we present in Fig.[B]DMRG results for 
the short-range correlations as a function of J3 (J4 = 0). 
The state at (J3, J4) = (1,0) is known to be gapped4£ 
The essential information in Fig. [5] is that the curves are 
devoid of any features that might suggest a change of the 
phase. Thus we can assert that the phase at J3 = 1 is 
smoothly connected to the phase at J3 = 0, which is a 
RS phase. The variation of the gap with J4 is shown 
in Fig. [7] We see that the gap goes to zero around the 
point J4 = 0.710 identified above as a phase transition 
point to another singlet phase. Below we discuss in more 
detail the structure of the low-lying excitations close to 
J4 = Jf . 



2. Haldane-like phase 

The discussed mapping of Eq. ([2]) on the frustrated 
spin-1/2 ladder model suggests that the HL phase should 
occupy some region in the phase diagram for J4 > JJ 1 . 
To reveal the peculiarities of the suggested HL phase - 
as compared to the well-known Haldane phase of the pe- 
riodic spin-1 Heisenberg chain - notice that in the sector 
[1, 1, . . . , 1] the Haldane state \^h) is the exact ground 
state of the mixed-spin Hamiltonian ([2]) at the symmetric 
point J3 = J4. In the general case (J3 7^ J4), the energy 
of this state Eh — hIH^ h) reads as 

where s H = -1.40148403897(4) is the the ground- 
state energy per bond of the periodic spin-1 Heisenberg 
chain4i Here, we have used the fact that the operator V 
[Eq. ([3])] does not have non-zero matrix elements in the 
sector [1,1,..., 1]: In particular, we have ("J/jy | V\^h) = 
0. The energy of the Haldane state Eh as a function of 
J4 ( J\ = Ji = J3 = 1) is shown in Fig. |4] (the ab line). 
Interestingly, at the special point J 4 = J| 2 = 0.875 - also 
related to an abrupt change of the spin-1 rung correla- 
tions - the DMRG estimate for the ground-state energy 
of the Hamiltonian @ E/L = —1.6899 almost coincides 
with the energy of the Haldane state (Ejj/L = —1.6889) 
obtained from Eq. (fTTj) . As already mentioned above, 
the numerical analysis implies that the special point J| 2 
is a quantum phase-transition point from a singlet non- 
degenerate state to a state exhibiting a net magnetic mo- 
ment. The above remarks suggest that the HL phase ap- 
pears as a good candidate for the phase diagram of the 
mixed-spin model. 

Further qualitative information about the characteris- 
tics of this phase can be extracted from a perturbative 
analysis starting from the symmetric point J3 = J4 and 
based on the Haldane state in a periodic spin-1 chain. 
Note that in some interval (J4 < J| 2 ) the parameter J a , 
which controls the V term in Eq. ((2|), may be used as 
a small parameter (e.g., J a = 0.0625 for J4 = 0.875). 
Thus, up to second order in J a , the ground-state energy 
takes the form E = Eh — const (1 — J 4) L, where const 
is some positive number of order one. Qualitatively, this 
result reproduces the behavior of the ground-state energy 
in the interval Jf 1 < J4 < Jf 2 extracted from the DMRG 
analysis (see Fig. [4}. To some extent, this result also val- 
idates the choice of \^h) as a starting unperturbed state. 

As compared to the Haldane state, some peculiarities 
of the HL phase can be revealed by looking at the first- 
order correction in J a to the wave function [^h), 

\9) = \*H)+JaY, { *E V fl£ ) ^ +0 W)- (12) 

Here the sum runs over the excited eigenstates \^ n ) of 
the Hamiltonian @ at J 3 = J 4 , and Eq = Eh- The 
matrix elements of V (see the Appendix) admit only two 
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types of excited states (^1,2) ) defined, respectively, in 
the sectors [f , . . . , f , 0, 0, 1, . . . , 1] (two neighboring rungs 
in singlet states) an [1, 1, . . . , 1, 2, 0, 1, . . . , 1] (one rung 
in a a quintet state an a neighboring rung in a singlet 
state). The weights of both types of defect configurations 
in the HL state change in the interval J| x < J4 < J| 2 : 
While the weight of the | M/i) configurations grows in a 
region around the transition point J% , the con- 
figurations (containing spin-2 defects) become visible in 
the DMRG result for the spin-1 rung correlations only 
in a short interval preceding the transition to a mag- 
netic state (see Fig. [3J . Note that the observed increase 
of the weight of the configurations formally con- 

tradicts the perturbation result in Eq. (fl~2"j) . which pre- 
dicts the opposite behavior. A reasonable resolution for 
this is provided by the guess that close to the transi- 
tion point J| 2 some of the eigenenergics E n related to 
the sector [1, 1, . . . , 1, 2, 0, 1 . . . , 1] soften. As of now we 
do not have firm numerical results in favor of such a sug- 
gestion, although some preliminary DMRG results, using 
open boundary conditions, seem to predict strong reduc- 
tions of the singlet-quintet and triplet-quintet gaps close 
to Jf. 
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FIG. 7: DMRG and ED numerical results for the singlet- 
singlet (A„) and singlet-triplet (A s ) excitation gaps in the 
mixed-spin model JlJ vs. J 4. The DMRG data points corre- 
spond to extrapolated values of A s obtained by a polynomial 
fit (up to L = 90) for open boundary conditions. The ED 
data concerns periodic clusters with L = 10, 12, 14. 



3. The RS-HL transition 



Turning to the region around the transition point J% 1 , 
it is instructive to comment on our numerical results 
for the excitation gaps (Fig. [7]) in the light of the dis- 
cussed mapping to the spin-1/2 ladder model. For the 
latter model, it has been numerically established 3 -^ that 
(i) the lowest state above the singlet ground states close 
to the phase boundary is a singlet excitation and (ii) 
the low-lying triplet excitations are gapped in the whole 
region of the phase diagram in Fig. [5j including the 
phase-transition boundary. Such a structure of the low- 
lying excitations is consistent with the established first- 
order quantum phase transition, which is described as 
a level crossing of two singlet ground states. As al- 
ready mentioned, the character of the RS-HL transition 
in the weak-coupling limit (Jj_, J 4 <C J 3 ) is still under 
debate^ ' 36 ' 37 As a matter of fact, there are some indica- 
tions for a second-order RS-HL transition^ and an inter- 
mediate dimer phas o 34 i 37 , but the debate concerns only 
the weak-coupling part of the phase boundary. Looking 
at the coordinates of the A and A images of the tran- 
sition point Jf 1 (Fig. [5]), it is clearly seen that the dis- 
cussed RS-HL transition at J 4 = J| x does not belong to 
the weak-coupling region. Hence, one may expect a first- 
order RS-HL transition at Jf related to a level crossing 
of singlet ground states. 

Figure [7] presents our numerical (DMRG and ED) re- 
sults for the singlet (A s ) and triplet (A t ) gaps of the 
lowest excited modes above both singlet ground states. 
Let us first discuss the ED data for the gaps. As clearly 
seen, both minima, related to the A s and A f data points, 
are located close to the expected transition point at 



J4 = 0.710. More importantly, an extrapolation of the 
ED data for J4 = 0.710 implies that the A s points scale 
to smaller values than A t . This observation is consistent 
with the expected low-energy structure close the first- 
order transition point between the RS and HL phases. 

Turning to the DMRG results for A t (J4), one observes 
that the triplet gap of the RS phase takes very small val- 
ues close to the suggested transition point (J4 = 0.710). 
We could not conclusively exclude the possibility of a 
gapless triplet excitation at the transition point. In any 
case, such a behavior indicates some peculiarities of the 
RS-HL transition in the mixed-spin system, as compared 
to the uniform-spin case. Another issue to be noticed 
is the steep (but definitely finite) slope of the function 
A t (J4) at the transition point. This suggests a relatively 
large correlation length of this triplet excitation close to 

TCl 

j 4 . 



C. Ferrimagnetic phases 

Looking at the DMRG results for the short-range cor- 
relations (Fig. [3]) , it is easy to realize that a ferrimagnetic 
phase, closely related to the ferrimagnetic ground state 
of an antifcrromagnctic Hcisenberg chain with alternat- 
ing (2, 1) spins, is stabilized around the symmetric point 
J4 = 1. Exactly at J 4 = 1, the ground state of the 
Hamiltonian ([2]) belongs to the sector [1, 2, . . . , 1, 2], so 
that both models are equivalent in the low-energy sector 
of the spectrum. The discussed ferrimagnetic phase (Fi) 
exhibits the magnetic moment per rung Mq = 1/2 and 
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(a) (b) 




Site index Site index 

FIG. 8: (a) The local magnetizations (sf^n) an d (ci,2n+i) 
(n = 1, . . . , 50) along the first leg as a function of the site 
index. The data shown is for Ja = 1.55. (b) The spin-1 
rung correlations along the length of the ladder (L = 100) at 
J 4, = 0.90. The values show a clear alternation between ps 1 
and « — 1 which indicates a two sublattice structure and a 
doubled unit cell containing four rungs. 



survives almost in the entire region after J| 2 , excluding 
some narrow interval in the vicinity of the latter point. 
This is also seen in Fig. |5Ja) which shows a typical be- 
havior of the local magnetizations (sf 2n ) and (ci2n+i) 
(n = 1, . . . , L/2) along the first leg at J4 = 1.55. The 
values of the spin-1 and spin- 1/2 magnetic moments are 
0.866950 and —0.366950, respectively. We see that the 
sum of the local magnetic moments is 1/2, as expected in 
a Lieb-Mattis type ferrimagnetic state with a quantized 
magnetic moment per rung Mq = 1/2. The deviations at 
the end are essentially because of open boundary condi- 
tions. We have verified numerically that these values do 
not change much after J4 = 1. 

For the region close to J| 2 , the DMRG results pre- 
sented in Fig. Elb) demonstrate the appearance of an- 
other ferrimagnetic phase (F2) in a narrow range of J4 
starting from the transition point J| 2 = 0.875 and ter- 
minating at J| 3 = 0.975. The F2 phase is characterized 
by the magnetic moment per rung Mq = 1/4. As clearly 
seen in Fig. [5Jb), in the F2 phase the space variation of 
the spin-1 rung correlations follow strictly the periodicity 
of the spin structure in the sector [2, 1, 1, 1, . . . , 2, 1, 1, 1]. 
Such a breaking of the translational symmetry is also seen 
in the inset of Fig. [31 where on the vertical axis we have 
plotted the magnitude of the difference of the spin-1 rung 
correlations in two neighboring unit cells for all values of 
J 4 . Clearly, the F 2 phase represents a two- fold degen- 
erate ground state, which is invariant under the transla- 
tion by two lattice periods. Our numerical analysis does 
not support the appearance of ferrimagnetic phases with 
larger periods. 



IV. CONCLUSION 

In conclusion, we have analyzed the combined effect 
of the quantum fluctuations and the competing inter- 
actions in a mixed-spin ladder composed of spin-1 and 
spin- 1/2 rungs which is closely related to a recently syn- 
thesized quasi- ID ferrimagnetic material. A comparison 



of the classical and quantum phase diagrams reveals the 
following changes in the related quantum system. As ex- 
pected, the classical ferrimagnetic phase also presents in 
the quantum phase diagram, but there appears another 
two-fold degenerate ferrimagnetic state which breaks the 
translational symmetry. As may be expected, the clas- 
sical Neel state does not survive quantum fluctuations. 
More interestingly, the classical canted state also com- 
pletely disappears. This is in contrast to some other ID 
spin systems exhibiting classical canted states where 
this type of classical magnetic order partially survives 
quantum fluctuations. In the present case, both the clas- 
sical long-range ordered states are replaced by two singlet 
non-degenerate gapped states (RS and HL). 

Turning to the weakly frustrated region, it has been 
established that the behavior of the system strongly re- 
sembles that of a two- leg spin- 1/2 Hciscnbcrg ladder with 
frustrating diagonal interactions. However, concerning 
the quantum phase transition between the RS and HL 
phases, we have found a few indications demonstrating 
some peculiarities (such as the extremely small triplet 
gap at the transition point) of the mixed-spin system. 
These issues deserve further investigations. 

Finally, although the available experimental results on 
the ferrimagnetic ladder material Fe 77 Fe 777 (trans-1,4- 
cyclohexanedicarboxylate) seam to point toward the re- 
alization of the Fi ferrimagnetic stated a detailed com- 
parison with the experiment requires a more extensive 
analysis of the quantum phase diagram including, e.g., 
different rung couplings J\ ^= J 2 , different pairs of rung 
spin magnitudes, and some anisotropics. Concerning the 
condition J3 = 1, as shown in Fig. [5] it simply restricts 
the path in the more general parameter space ( J3 ^ 1) to 
a straight line crossing one and the same phase bound- 
ary. Therefore, there should be a relatively large region 
with J3 7^ 1 showing the same structure of the phase 
diagram. As to the second restriction (J% = J 2 ), its re- 
moval may be generally expected to bring new quantum 
spin phases. However, in both cases we have numerically 
checked that relatively small deviations from the condi- 
tions Ji = J 2 = J3 do not bring qualitative changes on 
the established quantum phase diagram. 
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Appendix: Projection onto the spin- 1/2 ladder 

We have to project the spin-1 rung states onto the 
states of the spin-1/2 rungs. To this end, we use the 
projection operator P = P\Pi, . . . , Pl, where the rung 
projection operator P n reads as 

P n = J2\T? n )(T? n \, a = 0,x,y,z. (A.l) 

a 

Here \T 2n ) denotes the singlet state of the 2nth spin-1 
rung and |T|„) = (i / ' y/2)e Mm \l)\m) are the triplet states 
of the same rung in a vector basis which is a tensor prod- 
uct of the vector bases of the spin-1 objects (i.e., \x), \y), 
and \z)). In the following, the Greek indices take the 
values 0,x,y, and z, whereas the Latin ones - x,y, and 
z. 

Up to first order in J a , the projected Hamiltonian reads 

as 

n eff = PUP. (A.2) 



By using the expressions for the matrix elements 
ra4JT 2 « n ) = 25™, (Tl\L k 2n \Tl) = (TZ\L k 2n \Tl) = 
0, and (TZ\L k 2n \T° n ) = -2^2/35 mk , one obtains 

P n sj n P n =2j2\Tt)(TL\ = "In, (A.3) 

k 

where a%n is an effective rung-1/2 spin operator, and 
P V P — 

- 2 \/f E [\ T L)( T 2n\ + \tL)(tL\] (lL-l + lL + l) ■ 



Note that the operator in the square brackets is an ef- 
fective l2n rung vector operator for spin-1/2 rungs. Sum- 
ming the above results, we obtain the effective spin-1/2 
ladder model presented in Eq. (|9]). 
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